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ABSTRACT 

A significant and growing portion of systematic error on a number of fundamental parame- 
ters in astrophysics and cosmology is due to uncertainties from absolute photometric and flux 
standards. A path toward achieving major reduction in such uncertainties may be provided by 
satellite-mounted light sources, resulting in improvement in the ability to precisely characterize 
atmospheric extinction, and thus helping to usher in the coming generation of precision results in 
astronomy. Using a campaign of observations of the 532 nm pulsed laser aboard the CALIPSO 
satellite, collected using a portable network of cameras and photodiodes, we obtain initial mea- 
surements of atmospheric extinction, which can apparently be greatly improved by further data 
of this type. For a future satellite-mounted precision light source, a high-altitude balloon plat- 
form under development (together with colleagues) can provide testing as well as observational 
data for calibration of atmospheric uncertainties. 



Subject headings: techniques: photometric, instrumentation: photometers, atmospheric effects, dark energy, space 
vehicles: instruments, balloons 



1. INTRODUCTION 

While our understanding of the Universe has 
changed and improved dramatically over the past 
25 years, the improvement of our knowledge of 
absolute spectra and flux from standard calibra- 
tion sources, upon which the precision of mea- 
surements of the expa nsion history of the Universe 
(jAlbrecht et a/. || 2006l), and of s t ellar and galactic 
evolution [e.g. lEisenstein et all (|200ll )] are based, 
has been far slower and not kept up with reduc- 
tion of other major uncertainties. As a result, un- 
certainties on absolute standards now constitute 
one of the dominant systematics for measurements 
such as the expansion h istory of the Universe usin g 
type la supe r novae fe. o.lWood-Vasev et all ([20071 ). 



Astier et all (|2006l ). iKnop et all ([20031 )] . and a 



significant systematic for measuremen ts of stellar 



popu l ation in galaxy clu ster counts ([Kent et 



20091 : iKoester et aHhoOlh . and upcoming photo- 
metric redshift surveys measuring growth of struc- 



ture ( Connolly et al. 2006). There are prospects 
for improvement in uncertainties from standard 
star flux and spectra ([Kaiser et aZ.ll2008h . but the 
traditional techniques of measurement of standard 
stellar flux from above the atmosphere suffer from 
basic and inherant problems: the variability of 
all stellar sources, and the difficulty of creating 
a precisely calibrated, cross-checked, and stable 
platform for observation above the Earth's atmo- 
sphere. 

The presence of an absolute flux standard in 
orbit above the Earth's atmosphere could pro- 
vide important cross-checks and potential sig- 
nificant reduction of photometric and other at- 
mospheric uncertainties for measurements that 
depend on such calibration. Monochromatic 
sources, especially ones which could cover mul- 
tiple discrete wavelengths, or tune over a spec- 
trum, could potentially help to further reduce 
spectrophotometric error. For calibration of tele- 



1 



scope optics and de tector characteristics, authors 



(jStubbs et a/.ll2006l ) have both conceived of and 



used a wavelength-tunable laser within present 
and upcoming telescope domes as a color calibra- 
tion standard. Although a wavele ngth-tunable 
laser calibration source in orbit ([Albert et al. 
2006) does not exist yet, at present there is a 
532 nm laser in low-Earth orbit pointed toward 
the Earth's surface, with precise radiometric mea- 
surement of the energy of each of the 20.25 Hz 
laser pulses, o n the CALIPSO sat ellite, launched 
in April 2006 (| Winker et aLll2009h . We have col- 



lected data from a portable network of seven cam- 
eras and two calibrated photodiodes, taken during 
CALIPSO flyovers on clear days in various loca- 
tions in western North America. The cameras 
and photodiodes respectively capture images and 
pulses from the eye-visible green laser spot at 
the zenith during the moment of a flyover. Us- 
ing precise pulse-by-pulse radiometry data from 
the CALIPSO satellite, we compare the pulse en- 
ergy received on the ground with the pulse energy 
recorded by CALIPSO. The ratio determines the 
atmospheric extinction. 

These measurements apply only at the zenith, 
at the specific location of the telescope network, 
for the specific laser frequency, at the time of the 
CALIPSO flyover. To obtain results that are more 
easily applicable to modern astronomy, a satellite 
laser source must be able to point at major observ- 
ing facilities, and modify its frequency to be able 
to scan through the visible and near-infrared spec- 
trum. CALIPSO does not have this ability, how- 
ever future atmospheric science satellite missions 
may have such capabilities. A brief discussion of 
this is in Section 8. 

For the analysis of the images and the deter- 
minations of absolute flux, both laboratory and 
field calibration of the cameras and the photodi- 
odes was essential. We developed a calibration 
system to determine the sizes of effects such as 
the anisotropy, nonlinearity, and temperature de- 
pendence of the responses of the CCDs, so that we 
could accurately and precisely measure optical en- 
ergy of a source from its image. We were able to 
calibrate the cameras to a measured uncertainty 
of 2.5% or better in absolute flux. We include a 
detailed discussion of systematic effects and un- 
certainties on our measurements, currently domi- 
nated by atmospheric scintillation. Following the 



discussion of our results, we elucidate ways that 
this technique could be considerably improved be- 
yond this initial study. 

2. SATELLITE-MOUNTED LIGHT 
SOURCES 

Throughout history prior to 1957, the only 
sources of light above the Earth's atmosphere were 
natural in origin: stars, and reflected light from 
planets, moons, comets, etc. Natural sources 
have of course served extremely well in astronomy: 
through understanding the physical processes gov- 
erning stellar evolution, we are now able to pre- 
cisely understand th e spectra of st ars used as cali- 
bration sources [e.o. lBohlinl d2000h ]. Nevertheless, 
in all stars the vast bulk of material, and the ther- 
monuclear processes that themselves provide the 
light, lie beyond our sight below the surface of 
the star. Superb models of stellar structure are 
available, but uncertainties of many types always 
remain. 

Since the launching of the first man-made satel- 
lites, a separate class of potential light sources 
in space has become available. Observable light 
from most satellites is primarily due to direct so- 
lar reflection, or reflection from Earth's albedo. 
While providing a convenient method of observ- 
ing satellites, this light is typically unsuitable for 
use as a calibrated light source due to large un- 
certainties in the reflectivity (and, to a lesser ex- 
tent, the precise orientation and reflective area) of 
satellites' surfaces. Reflected solar light has, how- 
ever, been successfully used as an absolute infrared 
calibration source by the Midcourse Space Exper- 
iment (MSX), using 2 cm diameter black-coated 
spheres ejected from the MSX satellite, whose in- 
frared emissio n was monitored b y the instruments 
aboard MSX (jPrice et aZ.I l2004h . This technique 



proved highly effective for the MSX infrared cal- 
ibration; however, the technique is not easily ap- 
plicable to measuring extinction of visible light in 
the atmosphere. 

Many satellites have retroreflective cubes in- 
tended for use in satellite laser ranging. Reflected 
laser light from retroreflectors is critical for dis- 
tance measurements using precise timing; how- 
ever, like solar reflection, retroreflected laser light 
unfortunately also suffers from uncertainties in the 
reflectivity of the cubes, and in reflectivity as a 
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function of incident angle, that are too large to 
provi de a means of measuring atmospheric extinc- 
tion (jMinotdll974j ). A satellite-based mirror, ei- 
ther planar or spherical, which could reflect ei- 
ther sunlight or a ground-based light source, would 
also similarly require some means of calibrating 
reflectivity, both as a function of incidence an- 
gle and as a function of time (due to pitting from 
dust particles intersecting its orbit) — although it 
is a potentially interesting concept. Thus we are 
left with dedicated light sources aboard satellites 
themselves as the practical means of having a 
satellite-based visible light source for precise cali- 
bration of ground-based telescopes. 

Many satellites also carry some means of pro- 
ducing observable visible light, for self-calibration 
purposes or otherwise. The Hubble Space Tele- 
scope is one of many satellites carrying tung- 
sten, as well as deuterium , lamps for absolut e 



self-calibration purposes (|Pavlovskv et all [2001) 



Lamps for self-calibration are not limited to 
space telescopes for astronomy; earth observa- 
tion and weather satellites also commonly use 
i nternal tungsten lamps a s calibration sources 



( Nit hianandam et ali 119931 ). Such internal cali 



bration lamps are typically limited by the fact 
that they can degrade individually, and can be 
compared only with astronomical sources after 
launch, leaving stellar light as the only practical 
way to "calibrate the calibration device." Thus, 
such devices typically provide a cross-check rather 
than the basis for a true absolute irradiance cali- 
bration, or provide a means f or a separate calibra- 



tion, such as flat-field [e.g. iBohlin and Gillilandl 



(|2004l )]. Furthermore, present-day internal cali- 
bration lamps aboard satellites are certainly not 
intended for, nor are capable of, a direct cali- 
bration of the atmospheric extinction that affects 
ground-based telescopes. 

However, a satellite-based absolute calibration 
source for ground-based telescopes is not tech- 
nically prohibitive. As an example, a standard 
household 25-watt tungsten filament lightbulb 
(which typically have a temperature of the or- 
der of 3000 K and usually produce approximately 
1 watt of visible light between 390 and 780 nm) 
which radiates light equally in all directions from 
a 700 km low Earth orbit has an equivalent bright- 
ness to a 12.5-magnitude star (in the AB system, 
although for this approximate value the system 



makes little difference). In general, the apparent 
magnitude of an orbiting lamp at a typical incan- 
descent temperature which radiates isotropically 
is approximately given by 



-5.0 log- 



On (i^t)Y 
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5.9, 



(1) 



where P is the power of the lamp in watts, and h 
is the height of the orbit in kilometers. 

The dominant uncertainties in the amount of 
light received by a ground-based telescope from 
such an orbiting lamp would stem from uncertain- 
ties in onboard radiometry monitoring the power 
output of the lamp, any unsubtracted background 
from reflected earthshine, sunshine, moonshine, or 
starlight from the surface of the satellite itself, and 
potential deviations from perfect isotropic output 
(i.e. differences in the onboard-monitored output 
vs. the ground-observed output) of the light from 
the lamp. The uncertainty on the magnitude 
of the lamp stemming from uncertainties in the 
radiometrically-monitored output power would be 
limited by the precision of current radiometer 
technology. Modern radiometers, using electrical 
substitution radiometry, can achieve a precision 
of approximately 200 parts per million, with the 
domi nant uncertainty b eing the size of the aper- 
ture (jKopp et al. 120051) . The uncertainty on the 
magnitude of the lamp due to unsubtracted back- 
ground from reflected earthshine, sunshine, moon- 
shine, or starlight from the surface of the satellite 
would depend both on the size and surface ma- 
terial of the satellite, and on the performance of 
shuttering of the lamp to provide images to sub- 
tract such background. To minimize reflectance, 
one could coat the surface of the satellite with a 
black nonreflective coating, however this would re- 
sult in heating of the satellite as it passes through 
sunlight, until such heating was in equilibrium 
with thermal radiation from the satellite, an av- 
erage temperature of approximately 280 K for a 
spherical nonreflective satellite in continuous sun- 
light that generates no internal power, and 280 
K x (1 + 9.3 x 10~ 4 ( J^)) 1/4 for such a spherical 
nonreflective satellite of diameter d meters that 
generates P watts of internal power. The tem- 
perature differences as the satellite passes through 
sunlight and through the Earth's shadow would 
be quite large, and power supplies for a lamp 
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(and the lamp itself) would need to allow for such 
an operating temperature shift. Assuming this is 
achieved, there would be residual reflectance from 
the black coating. Typical NiP black coatings can 
achieve a reflectance of as l ow as 0.2% throughou t 



most of the visible range ( Kodama et al. 1990) 



Assuming perfect shuttering of the lamp output, 
the image-to-image variance in reflected light as 
the satellite passes through Earth's shadow would 
form the uncertainty, a value that is likely small 
and could be calibrated. Finally, the uncertainty 
on the magnitude of the lamp due to potential de- 
viations from perfect isotropic output of the light 
from the lamp would be able to be tightly con- 
strained using measurements done in the labora- 
tory before launch, as well as minimized by plac- 
ing the lamp inside an integrating sphere, and thus 
should be a very small, if not a negligible, contri- 
bution. Remaining uncertainties from this source 
could stem from uncertainty in the precise orien- 
tation of the satellite, and/or effects from degra- 
dation of the lamp on the isotropic output of the 
light, but would likely be small. Thus the pre- 
dicted total sysmatic uncertainty on the radiance 
of an optimally-designed orbiting lamp would be 
dominated by the precision of radiometric moni- 
toring technology, and be in the range of approxi- 
mately 200 parts per million if the best presently 
available radiometric technology is used. 

An alternative to an isotropic or near-isotropic 
lamp would be a laser source, with beam pointed 
at the observer (with a small moveable mirror, for 
example). Divergences of laser beams are typi- 
cally on the order of a milliradian (which can be 
reduced to microradians with a beam expander) 
so much less output power than a lamp would be 
required for a laser beam to mimic the brightness 
of a typical star. (However, note that the wall- 
plug power of typical diode-pumped lasers is typ- 
ically in th e vicinity of 20 tim es the output laser 
power [e.g. ISeas et all ( 2007 )]. which is further- 



more a large reduction in wall-plug power from 
the 10 4 level typical of flashlamp-pumped lasers.) 
The apparent magnitude of an orbiting laser with 
Gaussian beam divergence pointed directly at a 
ground-based telescope, is given by 



m 



-2-51og 10 



P 



20.1, 



(2) 



where P is the laser power in milliwatts, h is the 



height of the orbit in kilometers, and d is the RMS 
divergence of the laser beam in milliradians, un- 
der the assumption that the aperture of the tele- 
scope is small compared with the RMS width of 
the beam at the ground, hd. The RMS divergence 
would be the combination of the divergence at the 
source, and the divergence due to the atmosphere. 
In clear conditions, total atmospheric divergence 
in a vertical pa th is at the ley el of approximately 
5 microradians (|Tatarskilll96ll ). and this of course 
only acts on the last fraction of the laser path that 
is within the atmosphere, so as long as the source 
divergence is significantly larger than this, atmo- 
spheric divergence would be negligible. Clearly, 
even a sub- milliwatt laser in low Earth orbit would 
need to have either its divergence increased at the 
source, or its power reduced via filtering, for it to 
be suitable for astronomical calibration. We shall 
consider the filtering option. The major uncertain- 
ties in the amount of light received by a ground- 
based telecope from an orbiting laser would stem 
from uncertainties in the pointing and beam pro- 
file of the laser light (which would likely need to 
be monitored by an array of small dedicated tele- 
scopes outboard of the main ground-based tele- 
scope), in time-dependent variation of the laser 
output power (which would likely need to be moni- 
tored by onboard radiometry) , and in degradation 
of the filter over time. Nevertheless, laser light 
has the benefit of being monochromatic, allowing 
for calibration of individual wavelengths. With a 
widely-tunable laser, an entire spectrum could be 
calibrated, removing the significant inherant un- 
certainties associated with comparing the spectra 
of astrophysical objects with the spectrum of a 
calibration source. 

The uncertainty on the apparent magnitude of 
an orbiting laser stemming from uncertainties in 
the radiometrically-monitored laser power would 
be limited by the precision of current radiome- 
ter technology. Modern electrical substitution ra- 
diometers can achieve a precision of approximately 
100 parts per million when aperture uncertain- 
ties can b e neglected, as in the case of laser ra- 
diometry (|Kopp et aLll2005j ). Uncertainties on the 
magnitude due to uncertainty in the pointing and 
beam profile would potentially be limited by the 
size of the array of outboard telescopes for moni- 
toring the laser spot, and by calibration differences 
between the individual telecopes in the array and 
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with the main central telescope. The latter could 
clearly be minimized by a ground system for en- 
suring the relative calibration of the outboard tele- 
scopes and main telescope are all consistent. Thus 
one could achieve a similar, or potentially even 
smaller, uncertainty on apparent magnitude from 
an orbiting laser source as from an orbiting lamp 
source, but with the major added complication of 
the need for a large, dense array of outboard tele- 
scopes to very precisely monitor the position of 
the centroid of the laser beam on the ground. 

The uncertainties considered above assume that 
the exposure time is long compared with the co- 
herence time of the atmosphere. With short expo- 
sures, or in the case of a laser that either quickly 
sweeps past, or is pulsed, atmospheric scintilla- 
tion can play a major role in uncertainty in ap- 
parent magnitude of a satellite-mounted source. 
A typical timescale for a CW laser with 1 millira- 
dian divergence in low Earth orbit to sweep past 
is tens of milliseconds, which is of the same order 
as characteristic timescales of atmospheric scin- 
tillation, and the typical timescale of single laser 
pulses is nanoseconds, much shorter than scintilla- 
tion timescales, thus one cannot assume that such 
effects can be time-averaged over. In idealized con- 
ditions, for small apertures D <~ 5 cm and sub- 
millisecond integration times, the relative stan- 
dard deviation in intensity 07 = A.I/{I), where 
AI is the root-mean-square value of /, is given by 
the square root of 



erf = 19.12A~ 7/6 / 
Jo 



C*(h)h 5/e dh, (3) 



where A is optical wavelength (in meters), C^(h) 
is known as the refractive-index structure coeff i- 
cient, and h is altitude (in meters) ( Tatarskll96ll ). 
Large apertures D >~ 50 cm have a relative stan- 
dard deviation in intensity given by the square 
root of 



f c 

aj = 29.48L>- 7/3 / 
Jo 



Cl{h)h 2 dh (4) 



(|Tatarskilll96lh . The values and functional form 
of Cn(h) are entirely dependent on the particular 
atmospheric conditions at the time of observation, 
however a relatively typical profile is given by the 
Hufnagel- Valley form: 



Cl{h) = 5.94 x 10- 53 (t./27) 2 / l 10 e-' 1 / 1000 

-16 -h/lbOO 



wher e A and v are free parameters ([Hufnagel 
1974). Commonly- used values for the A and v pa- 
rameters, which represent the strength of turbu- 
lence near ground level and the high-altitude wind 
speed respec tively, are A = 1.7 x 10~ 14 m~ 2 / 3 and 
v = 21 m/s ( Roggemann and Welsh 19961 ). Using 
these particular values, for a small aperture, the 
relative standard deviation <jj would be expected 
to be 0.466 for 532 nm light, which is not far off ex- 
perimental scintillatio n values for a c l ear ni ght at 
a typical location [e.g. Jakeman et all ( 1978j )]. For 



a single small camera, this is an extremely large 
uncertainty. Other than by increasing integration 
time (which is not possible with a pulsed laser) 
or by significantly increasing the camera aper- 
ture, the only way to reduce this uncertainty is 
to increase the number of cameras. With N cam- 
eras performing an observation, which are spaced 
further apart than the coherence length of atmo- 
spheric turbulence (typically 5 to 50 cm), the un- 
certainty from scintillation can be reduced by a 
factor y/N (for large N). 

The analysis above considers a hypothetical 
pointable satellite-mounted calibration laser, and 
is necessarily both speculative and approximate. 
However, at present there is an actual laser in 
low Earth orbit, visible with both equipment and 
with the naked eye, and analysis of ground-based 
observational data of the laser spot can be used 
for comparisons with the above, as well as for de- 
velopment of and predictions for potential future 
satellite-based photometric calibration sources of 
ground telescopes. 

3. THE CALIPSO SATELLITE AND 
GROUND-BASED OBSERVATION 
NETWORK (2006-08) 

The CALIPSO (Cloud Aerosol Lidar and In- 
frared Pathfinder Satellite Observations) satel- 
lite was launched on April 28, 2006 as a join t 
NASA and CNES mission (j Winker et al\\200$ i. 
CALIPSO is part of a train of seven satellites 
(five of which are orbiting at the date of this arti- 
cle), known as the "A- Train," in sun-synchronous 
orbit at a mean altitude o f approximately 690 
~ 120081 ). A typical 1-day 



km (Savtche nko et al. 



2.7 x 10" 



At 



-fc/100 



(5) 



ground track for CALIPSO is shown in Fig. [T] 
CALIPSO completes an orbit every 98.4 min- 
utes (approximately 14.6 orbits per day), and 



5 




Fig. 1. — A typical 1-day ground track of the CALIPSO satellite (|Garber et adl2007t) 



repeats its track every 16 days. CALIPSO con- 
tains a LIDAR (Light Detection and Ranging) 
system, known as CALIOP (Cloud Aerosol Li- 
dar with Orthogonal Polarization), with a pri- 
mary mission of obtaining high resolution vertical 
profiles of clouds and aeros ols in the Earth's at- 
mosphere (|Hunt et a/.ll2009h . The CALIOP laser 
produces simultaneous, co-aligned 20 ns pulses of 
532 nm and 1064 nm light, pointed a small angle 
(0.3°) away from the geodetic nadir in the forward 
along-track direction, at a repetition rate of 20.16 
Hz. The light enters a beam expander, following 
which the divergence of each laser beam wave- 
length is approximately 100 /irad, producing a 
Gaussian spot of approximately 70 m RMS diam- 
eter on the ground. The pulse energy is monitored 
onboard the satellite, and averages approximately 
110 mJ, at each one of the two wavelengths, per 
pulse. The effective apparent magnitude of the 
532 nm laser spot at the precise center of the 
beam is thus approximately -19.2, however this 
high brightness, of course, falls off rapidly as one 
moves away from the center of the beam. 

During 2007, the CALIPSO beam was observed 
at several locations in western North America us- 
ing a portable ground station consisting of seven 
simultaneously-shuttered Panasonic FZ50 digital 



cameras, connected as shown in Fig. [3J and two 
calibrated photodiodes co-located with the central 
camera. One of the two photodiodes is a Hama- 
matsu S2281 silicon photodiode supplied by the 
National Institute of Standards and Technology 
(NIST, in Gaithersburg, MD), and calibrated by 
NIST to an absolute radiometric standard in the 
visible and near-infrared. A 12.5 mm diameter 
Semrock 532 nm line filter was mounted over the 
front face of this photodiode. The other photodi- 
ode is a Hamamatsu S1336-5BQ, with no filter. 

A list of the observations can be found in Ta- 
ble [TJ Two example observations can be seen in 
Figs. S] and [5] Observing locations were selected 
by means of the CALIPSO ground track, mon- 
itored by NASA. For each observation, a point 
along the ground track was selected the day prior 
to the night of the overpass by means of the follow- 
ing criteria: road accessibility, lack of obscuring 
trees and streetlights, avoidance of fenced-off pri- 
vate land, and local weather conditions. A desert 
environment is ideal; for this reason, several of 
the observations were performed in the southwest- 
ern U.S. Even thin cloud cover can create a differ- 
ence in time-averaged optical density between the 
CALIPSO beam path and light from stellar stan- 
dard sources on the image nearby; thus all obser- 
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Fig. 2. — Initial single-camera ground observation of a CALIPSO laser pulse, taken at 07:08:30 UTC on 
Nov. 23, 2006 near Granby, Quebec. 



vations were performed on clear nights. 

Each camera was placed pointing vertically, 
with the images approximately centered on the 
zenith, for all observations. Tripods were used for 
all observations except for the first (Apr. 17), in 
order to reduce issues of dust, or condensation, 
obscuring the camera lenses (for the Apr. 17 ob- 
servation, the cameras were placed directly on the 
ground). The exposure time for all observations 
with all cameras is 60 s., allowing sufficient expo- 
sure to capture images of nearby stars as well as 
the CALIPSO beam. For each camera, the aper- 
ture f-stops were set at f/4.0. As the lenses are 55 
mm focal length, the aperture diameters were thus 
13.75 mm. The images sensors for Panasonic FZ50 
cameras are 1/1.8" optical format CCDs (5.319 



mm x 7.176 mm) with 10.1 million effective pixels. 
The field of view of each camera is 20.5 x 15.5 de- 
grees; the axes are randomly oriented with respect 
to RA and declination. The two photodiodes, co- 
located with the central camera, were also pointed 
toward the zenith and each surrounded by a black- 
ened 1 cm radius tube, extending 10 cm above the 
face of the photodiode, in order to reduce back- 
ground light. Approximately 90 minutes prior to 
each overpass, a test image was taken with all 
seven cameras plus the two photodiodes, to ensure 
all observing elements were properly functioning. 
CALIPSO beam exposures were started approxi- 
mately 30 seconds before the overpasses. About 
half of all observation attempts were unsuccessful 
in obtaining an image of the CALIPSO beam in 
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Fig. 3. — CALIPSO satellite observation setup using seven Panasonic FZ50 cameras connected by coaxial 
cable, to obtain images of four pulses of the CALIPSO laser, and obtain constraints on the energy, position, 
and width of the beam pulse(s) on the ground. 



all cameras. The most common reason for failed 
attempts was uncertainty in the CALIPSO ground 
track. The uncertainty for CALIPSO ground track 
is approximately 100 meters, however this is the 
uncertainty of the a postiori CALIPSO ground 
track data (available three days after it is taken); 
for the a priori ground track forecast, the uncer- 
tainty is significantly greater, and is often subject 
to systematic biases. Thus there is a very large 
chance of missing the beam. In addition, in the 
few days following drag makeup maneuvers by the 
CALIPSO satellite (which are performed approxi- 
mately once per month) , the ground track forecast 
is especially subject to large systematic offsets. 
The second most common cause of failed observa- 
tions was condensation, or frost, on camera lenses; 
this was also a difficult problem to ameliorate for 
cameras with lenses that are neither shielded nor 
heated, other than by only choosing sites in which 
humidity is very low. 

4. CAMERA AND PHOTODIODE DATA 
FROM 2006-08 

The data taken by the Panasonic cameras is 
recorded in a proprietary Panasonic format (re- 



ferred to as a raw image). After transferring the 
data to computer disk, the raw images were con- 
verted into 16-bit TIFF images (using the sRGB 
colorspace) using the Silkypix Developer Studio 
2.0 SE program (v. 2.0.14.10). Then, using the 
GIMP 2.0 program, the TIFF images were sepa- 
rated into their red, green, and blue components, 
and each component color in each image was con- 
verted to a FITS image file. Using the DS9 pro- 
gram (v. 5.51), the approximate pixel location of 
the centroid of the laser spot on each image con- 
taining the green pixel data was determined man- 
ually. Using SuperMongo (v. 2.4.34), the pixel 
values in a 51 pixel by 51 pixel square centered 
around the approximate centroid pixel location 
were written out to a text file. 

The laser spots are easily found on the original 
color images by eye-scanning for a green-colored 
dot that has no trail. (The minute-long exposures 
are enough to create visible trails on stars and 
other celestial objects, whereas the laser flash, be- 
ing essentially instantaneous, creates no trail.) In 
all cases where a laser spot existed, it was obvious 
by eye, and there never was a significant question 
of whether it was the laser, or a green-colored ce- 
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April 17 observation: near Carefree, AZ, USA, at 09:34:30 UTC 




Fig. 4. — Seven-camera observation of a CALIPSO overpass, taken on Apr. 17, 2007 near Carefree, Arizona. 



lestial object. In the FITS files corresponding to 
the red and blue parts of the images, little or no 
signal from the CALIPSO spots is found, and thus 
the red and blue components add very little infor- 
mation to the measurement of the signal from the 
CALIPSO beam. Thus only the green parts of the 
images are used in the final analysis. 

The data from the two Hamamatsu photodi- 
odes was amplified by means of a Hamamatsu 
C9052 Si photodiode readout board in the light 
integration circuit configuration. The output was 
then subtracted for constant background light and 
dark current contribution, the former of which var- 
ied accoring to the amount of moonlight, etc., in 
the background, and thus was adjustible on-site 
by means of a potentiometer. The background- 



subtracted signal was then further amplified by a 
factor of 50. The signal was then digitized and 
read out at 200 Hz over the 60 s interval to com- 
puter disk by means of a LoggerPro™ readout 
device; then converted to a text file. Power was 
supplied via an AC inverter and automobile en- 
gine. 

Only 1 of the 6 camera observations (the one 
taken on 30 May) has properly-recorded photo- 
diode data. Photodiode readout electronics, or 
other, problems prevented successful recording in 
the other 5 observations. While the photodi- 
ode signal was visible in the 30 May observation, 
we observed significant differences in the ratio of 
photodiode-to-camera signal size from that which 
we measured in the laboratory, which could be due 
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May 1 Observation: west of the Great Salt Lake, UT, 9:44:50 UTC 




Fig. 5. — Seven-camera observation of a CALIPSO overpass, taken on May 1, 2007 near the Great Salt 
Lake, Utah. 



either to background moonlight affecting the pho- 
todiode signal, or differences in the gain of the 
photodiode electronics from that measured in the 
laboratory due to the lower temperature. Thus, 
we chose not to use the photodiode data which 
was collected in the field, even that from the 30 
May observation, in the analysis. Nevertheless, as 
the cameras were calibrated to the NIST photodi- 
ode absolute standard in the laboratory, the data 
from the seven-camera network allows the mea- 
surement of the CALIPSO pulse energies reaching 
the ground. 



5. ANALYSIS 

5.1. Absolute Photometric Calibration of 
the Panasonic Cameras 

The Panasonic cameras were calibrated to an 
absolute radiometric standard through the use of 
the NIST-calibratcd Hamamatsu photodiode (see 
above) and a Pavillion Integration W532-10FS 
low-noise 532 nm laser with Newport neutral den- 
sity filters mounted in the beam to reduce the in- 
tensity by a factor of approximately 2 x 10 5 . The 
cameras were mounted in the intensity-reduced 
beam, inside a light-tight enclosure, and each ex- 
posed for 1 minute, as during CALIPSO observa- 
tions. The cameras were then swapped with the 
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Fig. 6. — Measurement of the response linear- 
ity of the central camera CCD. The response is 
linear to within 2.4% over the energy scale of in- 
terest (i.e. below the plateau due to pixel satu- 
ration). Note that the laser power is attenuated, 
by a constant factor of O(10 5 ), using neutral den- 
sity filters. The response of the other cameras has 
similar linearity. 



calibrated photodiode and readings were taken. 
The analyzed images were then compared with the 
photodiode data to obtain a radiometric compar- 
ison point. 

Measurements of the response linearity, an- 
isotropy, and ambient temperature dependence of 
the cameras were also performed. Linearity was 
measured by modifying the power of the laser, us- 
ing input power control hardware supplied by the 
laser manufacturer. The radiometric calibration 
process above was repeated at 10 different laser 
power settings for two of the cameras. Linearity 
was found to be maintained to ±2.4%, as shown 
in Fig. [6] The anisotropy of the response of the 
CCDs was measured by obtaining images as above 
with the laser spot at 900 different places on each 
of the cameras' CCDs. The anisotropy, taken to 
be the standard devation of those 900 measure- 
ments, was found to be 0.5%, as shown in Fig. [7] 
Temperature calibration was performed using one 
of the cameras inside a sealed refrigerator with 
a small hole in it to allow entrance of laser light 
(from inside a light-tight tube). Over a temper- 
ature range from 20° C down to 0°C, the image 
intensity was found to increase by the surprisingly 



Fig. 7. — Anisotropy measurement of the CCD 
of the west camera. The horizontal axes show the 
x and y pixel position; the vertical axis is flux, 
in arbitrary units. The deviations from a flat dis- 
tribution on the vertical axis are magnified by a 
factor of 100. The other cameras have a similar 
degree of anisotropy (~0.5%). 



large value of 52%, as shown in Fig. [8] This vari- 
ation, alt hough large, is not an unexpe cted phe- 
nomenon (jHealev and Kondepudvlll994l ). 
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Fig. 8. — Measurement of the temperature depen- 
dence of the observed flux from a stable, calibrated 
laser source in the central camera (other cameras 
are similar). 
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5.2. Measurement of Camera Positions 

Knowing the precise relative positions of the 
cameras is necessary for extrapolating the cam- 
era energy measurements into a measurement of 
the total laser pulse energy at the ground. Dur- 
ing each observation, the location of each camera 
was measured using GPS, as well as a surveyor's 
tape measure to determine the relative positions 
via triangulation. An example set of locations of 
the cameras, and their uncertainties, is given in 
Fig. m 

The camera positions are overconstrained by 
the combination of GPS and measured distance 
information. The measured distances in gen- 
eral provide tighter constraints than GPS, how- 
ever the GPS measurements are needed to set the 
centroid position, and overall rotation angle, of 
the network. To find the best-fit positions and 
their uncertainties, a x 2 minimization fit is per- 
forme d using the ROOT graphica l analysis pro- 
gram ( Brun and Rademakersl 1997 ). 



5.3. CALIPSO Laser Pulse Energies 

CALIPSO measures the energy of each of its in- 
dividual laser pulses. The pulse energy monitoring 
on CALIPSO consists of NIST-calibrated photo- 
diodes mounted on an integrating sphere, with 
pulse-by-pulse energy measurement with design 
absolute precision of ±2% over full o rbit, and rela- 
tive precision of better than ±0.4% ( Winker et al. 



2009). The pulse energies are recorded in the 
CALIPSO datasets available from the Atmo- 
spheric Science Data Center (ASDC, located at 
NASA LaRCfl 

which may be compared with 
the pulse energies measured at the ground to 
obtain measurements of atmospheric extinction. 
The CALIPSO satellite-measured pulse energies 
which correspond to each of the six ground station 
observations are given in Table [TJ 

5.4. Extinction Results 

The image data text file produced for each ob- 
servation from each camera, as described above, 
must be fitted to a function that describes the 
CALIPSO laser spot as well as background from 
scattered light etc. (only rarely did any stars hap- 
pen to come within the 51 pixel square, corre- 



1 http:/ /eosweb. larc.nasa.gov 



sponding to about 0.1 square degrees, that con- 
tained the CALIPSO spot). The data were found 
to be consistent with a single 2-D Gaussian de- 
scribing the spot, plus a uniform background. 
However, for images in which the spot is especially 
bright, the center pixels often reach the satura- 
tion limit of the camera. Thus, to fit the data, a 
likelihood function that allows for pixel saturation 
(merely constraining the fitted function to be any- 
thing above the saturated pixel values, rather than 
the values themselves, in the case of the pixels that 
have reached the saturation limit — and other- 
wise being a standard % 2 likelihood) must be used. 
When combined with the Gaussian plus uniform 
background, this likelihood function describes the 
data well. 

For each image data text file, a total of five vari- 
ables ar e floated in the fit, which is p erformed in 
ROOT (jBrun and Rademakersl Il997l i: the x and 
y positions of the spot centroid, the standard de- 
viation of the Gaussian spot, the level of uniform 
background, and the total amount of light in the 
CALIPSO spot. Asymmetric Gaussian uncertain- 
ties are determined for all variables. 

Plots of pixel data in the region surrounding 
the CALIPSO spot in eight of the camera images 
can be found in Fig. [TO] Further plots and data 
may be found online at . 

The above fits provide a measurement of the 
amount of light recorded in each camera, in the 
digital pixel units of the camera. To translate 
these digital unit measurements into SI units of 
energy, we use the absolute photometric calibra- 
tion data for each camera as described above. This 
provides us with absolute energy measurements for 
each camera (with additional uncertainty due to 
the calibration- related error budget). 

The above procedure determines the amount of 
light reaching each individual camera. To deter- 
mine the total amount of CALIPSO light reach- 
ing Earth in a given observation, the data from all 
seven cameras must be combined together (along 
with the fitted positions of the cameras as de- 
scribed above). Thus, a secondary fit, again using 
ROOT, is performed which combines this infor- 
mation into a measurement of the total amount 
of light from a CALIPSO pulse reaching Earth. 
In this fit, the fact that more than one pulse 
may contribute to the light measured in the seven 
cameras must be taken into account. A func- 
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Camera Locations for 30 May 2007 CALIPSO 
Observation near Woss, BC, Canada 



Note: Satellite imagery was 
taker in 2003. At some point 
between then and 2007, the 
area was logged, and was 
free of trees when the 
observation took place. 
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Fig. 9. — Camera locations for the 30 May 2007 observation (as an example). Lengths were measured with 
a surveyor's tape measure, and uncertainties are estimated. Elevations are relative to the central camera 
and are estimated by eye (their uncertainties do not contribute greatly to uncertainties on the horizontally 
projected positions). Uncertainties on the GPS positions are approximately 3 m (the measured lengths 
provide tighter constraints). Between 2003 and 2007, the area was logged, and a small logging road had 
been cut along the locations of the east, west, and central cameras. 
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Table 1 

List of portable ground station observations of the CALIPSO beam in 2007 (with their 

GEOGRAPHIC POSITIONS AND TEMPERATURES AT THE TIMES OF OBSERVATION), THE RESPECTIVE 
MEASURED LASER PULSE ENERGIES FROM THE GROUND CAMERA NETWORK (-Eground) AND FROM THE 
CALIPSO SATELLITE ITSELF (-EcALIPSo), AS WELL AS THE MEASURED ATMOSPHERIC EXTINCTION 

RATIO (r e = Aground /-^CALIPSO)- 



Date Lat.t Long.* T obs (°C) M Aground (mJ) ^calipso (mJ) r e Notes 



13 Mar. 


49.090° N 


123.917° W 


5 ± 1 




111.3 ± 2.2 




(1) 


17 Apr. 


33.811° N 


111.996° W 


10 ± 3 


41.1 ± 24.1 


111.2 ± 2.2 


0.4 ± 0.2 


(2) 


28 Apr. 


37.606° N 


106.249° W 


-1± 3 




111.1 ± 2.2 




(3) 


1 May 


41.012° N 


112.917° W 


10 ± 5 


116.9 ± 52.4 


111.2 ± 2.2 


1 1+ - 4 


(4) 


30 May 


50.224° N 


126.561° W 


5 ± 3 


77.7 ± 37.4 


111.3 ± 2.2 


0.7 ± 0.3 


(5) 


1 Jun. 


49.123° N 


123.931° W 


11 ± 3 


108.8 ± 46.1 


111.4 ± 2.2 


1.0 ± 0.4 


(6) 



Notes: 

* Location of the central camera during the observation. 



^ T emperature data is from ( Environment Canadall2010l ) for the Canadian observations and from ( NOAA 
2010) for the U.S. observations. Uncertainties are author-estimated due to the fact that the weather stations 



are not located in the same places as the observations (e.g. the 1 May observation was nearly 200 miles from 
the nearest one) , nor do all the closest stations have data at the times of each observation. 

(1) Near Nanaimo, B.C., Canada. Seen only in 1 camera (thus no ground energy measurement) 

(2) Near Carefree, AZ, USA. Full 7 camera observation. 

(3) Near Monte Vista, CO, USA. Icing on camera lenses (thus no ground energy measurement). 

(4) Near the Great Salt Lake, UT, USA. Full 7 camera observation. 

(5) Near Woss, B.C., Canada. Full 7 camera observation. 

(6) Near Nanaimo, B.C., Canada. Full 7 camera observation. 



Table 2 

Systematic uncertainties (in % and in mJ) on the measured pulse energies. 



Systematic Uncertainty 


17 Apr. 


1 May 


30 May 


1 Jun. 




(% / mJ) 


(% / mJ) 


(% / mJ) 


(% / mJ) 


Atmospheric Scintillation^) 


46.7 / 19.2 


27.5 / 32.2 


32.7 / 25.4 


23.3 / 25.4 


CALIPSO Beam Profile Uncertainty^ 


35.0 / 14.4 


35.0 / 40.9 


35.0 / 27.2 


35.0 / 38.1 


Camera Throughput & Absolute Calibration^) 


5.0 / 2.1 


5.0 / 5.8 


5.0 / 2.9 


5.0 / 5.4 


Camera/ADU Statistics (a) 


1.7 / 0.7 


0.8 / 0.9 


2.4 / 1.9 


0.6 / 0.7 


Total 


58.6 / 24.1 


44.8 / 52.4 


48.2 / 37.4 


42.3 / 46.1 



Notes: 

As discussed in Section 6, the above uncertainties can be drastically reduced (to approximately 1.1%) by (a) 
the addition of more cameras, with larger aperture, (b) pre-flight precise measurement of beam profile and 
onboard monitoring, and (c) improvements in the laboratory calibration of the camera optics. 
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Fig. 10. — Camera image data for pixels in a 35 x 35 pixel region around recorded CALIPSO laser spots for 
eight example images. Pixel intensities are in camera analog-to-digital units (ADUs). The contours show the 
fitted 2-D Gaussian function. Pixels with values equal to the pixel saturation level (196) are considered in 
the fit to have true value either greater than or equal to this value (such pixels have \ 2 n °t constrained from 
above), as described in the text, thus the peak of the fitted function can often exceed the saturation value. 
There are 5 contours for each plot, equally spaced between the fitted function's maximum and minimum 
(usually 0) value in each histogram. 
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Fig. 11. — Camera time-integrated irradiance data, and the resulting fitted time-integrated irradiance 
maps, for the Apr. 17 (top) and May 1 (bottom) observations. The numbers at each camera refer to the 
measured value of time- integrated irradiance by that camera, with its associated uncertainty (68% CL), and 
the expectation value from the fitted function at the location of that camera. The contours on each plot 
are spaced at 1 /iJ/m 2 intervals. The upper-right inset on each plot extends the x and y axes in order to 
see the three 2-D Gaussians that comprise the fitted function (as described in the text), and the lower- left 
inset shows a different 3-D view (from the side, rather than from above) of the fitted function and the data 
points. 
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Fig. 12. — Camera time-integrated irradiance data, and the resulting fitted time-integrated irradiance maps, 
for the May 30 (top) and Jun. 1 (bottom) observations. Please see the caption of the figure on the previous 
page for a description of these plots. In both of these two observations, the location (foliage, etc.) prevented a 
second camera from being placed further east than the east camera, thus a "west-west" camera location was 
used instead of an "east-east" one. (Also note that central value of the measured time-integrated irradiance 
at the west- west camera location on the Jun. 1 observation is above the color scale chosen for that plot.) 
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tion consisting of three Gaussian pulses is fitted 
to the data from the seven cameras. (Note that a 
train of three pulses extends far further than the 
size of the camera network.) The direction head- 
ing of the CALIPSO satellite, and the separation 
between CALIPSO pulses, are known well from 
NASA ground track information and the precise 
20.25 Hz frequency of the laser respectively. How- 
ever, the centroid on Earth of a given pulse (in 
both latitude and longitude), as well as its total 
amount of light, must be fitted. Thus, those three 
variables are floated in the fit, and central values 
and asymmetric uncertainties are determined for 
each. The results of these fits are shown in Figs.fTTl 
and[H 

The measured energies of the CALIPSO pulses 
reaching Earth, denoted as -Eground, are given in 
Table □ The CALIPSO satellite-measured pulse 
energies, -EcalipsOj are also given, as well as their 
ratios r e which provide measurements of extinc- 
tion. 

Although the uncertainties are large, the re- 
sults, with the possible exception of the Apr. 17 
observation, are consistent with the expectation 
of r e as 0.95-1.0 (given the fact that all observa- 
tions were made on relatively clear nights). The 
low value for the Apr. 17 observation may be due 
to the fact that this observation was taken in a 
desert environment, and blowing dust might have 
partially obscured the camera lenses — especially 
since the cameras were not placed on tripods for 
that initial observation. 

The breakdown of the respective uncertain- 
ties on the measurements of energy reaching the 
ground are given in Table [2] As one can see, 



atmos pheric scintillation, calculated as per iFried 
(|1967l ). and uncertainty in the profile shape of the 
CALIPSO beam [the beam is assumed to be a 2- 
D Gaussian in shape; the uncertainty on the true 
shape is taken to represent a ±35% uncertainty o n 



the measured energy value ( Winker et a l. 2009)] 



are by far the dominant uncertainties. The for- 
mer can be reduced by increasing the aperture 
and the number of cameras, and precise pre-flight 
measurement of beam profile (along with onboard 
monitoring) can dramatically reduce the latter 
uncertainty. Increasing the number of cameras 
(and/or the bit depth of their CCDs) can also re- 
duce the small uncertainty due to ADU statistics. 
The other relatively small uncertainty due to cam- 



era throughput and absolute calibration (taken to 
be ±5% of the measured energy value) can be re- 
duced by functionality and use of calibrated photo- 
diode measurements during observations, as well 
as with more thorough laboratory studies of the 
cameras. 

One could compare these results with indepen- 
dent extinction measurements based on the mea- 
sured light from stellar standard sources within 
the camera images during observation, to de- 
termine their consistency. However, in light of 
the very large atmospheric scintillation and beam 
profile-related uncertainties on the CALIPSO- 
based measurements, such a comparison would be 
more fruitful with a future more precise and larger 
telescope observation network, and the next- 
generation LIDAR satellite following CALIPSO. 

6. MONTE CARLO STUDIES OF OPTI- 
MAL TELESCOPE NETWORK DE- 
SIGN 

In order to reduce the dominant uncertainty 
due to atmospheric scintillation, significantly 
larger apertures than provided by the small Pana- 
sonic cameras are required for aperture averaging 
over atmospheric cells. In 2009-10, a study was 
done to optimize the design of an entirely new 
network of telescopes and cameras to provide the 
best possible measurements of CALIPSO laser 
light energy reaching the Earth, at a given equip- 
ment cost. A Monte Carlo simulation program to 
perform this optimization was developed. 

Consumer 16" aperture Dobsonian reflector 
telescopes are available at relatively low cost, 
and can provide an accessible means of aperture- 
averaging over atmospheric scintillation with a 
pulsed laser satellite source. Compared with a 
relative standard deviation of observed energy 
due to atmospheric scintillation of 07 = 0.466 for 
small aperture cameras (such as the Panasonics) 
with pulsed 532 nm light (as described in Sec. 2), 
the relative standard deviation for a 16" ap erture 
telesco pe is 07 = 0.035 [calculated as per IFried 
(|1967l )]. a reduction factor far outweiging the in- 
crease in cost. 

We consider the optimization, for a given obser- 
vation attempt of a CALIPSO pulse, of telescope 
locations in a network of 20 such telescopes. The 
uncertainty on the predicted path of CALIPSO is 
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Fig. 13. — Optimized locations for the network of twenty 16" telescopes. The CALIPSO pulse locations are, 
of course, only expectation values: the location of the pulses has a Gaussian uncertainty of 50 m transverse to 
the satellite direction of motion, and the phase of the pulses in the train (the longitudinal location) is entirely 
uncertain (only the distance between the pulses is known) . The optimum length of the telescope network (as 
shown in the optimization plots in Fig.lMf is less than the separation distance between CALIPSO pulses due 
to the importance of having as many telescopes as possible observing a pulse (despite the relatively small 
risk of having the entire network fall in the gap between pulses in a given observation) . 



approximately 50 meters in the direction orthog- 
onal to the satellite's travel (which, for the pur- 
poses of this and the next paragraph, we will de- 
note by the y direction, with the direction of the 
satellite's travel denoted as the x direction). We 
take that uncertainty to be Gaussian distributed. 
CALIPSO pulses are separated by 333 m. In the x 
direction, the location of the pulses is entirely un- 
certain a priori; only their separation is knownj^l 
Thus no advantage is gained by having an ellipti- 
cal distribution of the telescopes, centered around 
some point in the x direction, rather than a rect- 



2 In principle, this may be able to be predicted to some de- 
gree by studying the timing of the pulses from 3 days earlier 
(CALIPSO data takes 3 days to become available for anal- 
ysis) and extrapolating via the 20.25 Hz pulse rate, but it 
is likely that this would not provide precise information on 
the position of the pulses. 



angular distribution. Thus, we consider a grid of 
telescopes, either a 4 x 5 rectangular pattern or 
4x5 hexagonal pattern {e.g. such as that shown 
in Fig. n~5|) . with 4 telescopes in the y direction. 
(Other arrangements, such as 5 x 4 or 2 x 10, were 
also considered; these are considerably less effec- 
tive.) In the y direction, the optimal placement 
of the telescopes is such that they are evenly dis- 
tributed in the space corresponding to the integral 
of the probability density distribution of the irra- 
diance of the beam pulses. The probability density 
of the satellite path is a Gaussian distribution with 
a = 50 m, and the pulses themselves are Gaussian 
with a = 35.25 m, thus the probability density 
distribution of the y projection of the irradiance 
is also Gaussian with a equal to the addition in 
quadrature of 50 and 35.25/V2, i-e. 55.87 m (with 
the factor of y/2 due to the projection of the 2-D 
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Gaussian pulses onto the 1-D y axis). The opti- 
mal placement of the telescopes will be where the 
integral of that Gaussian reaches 1/8, 3/8, 5/8, 
and 7/8; thus creating an even spread over the 
normal distribution. This occurs when the tele- 
scopes are 17.8 m and 64.3 m to the right and left 
(in the y direction) of the expectation value of the 
satellite path, corresponding to 46.5 m and 35.6 m 
separation between the telescopes in y, as shown 
for example in Fig. Q2J This sets the y distribu- 
tion. However, the placement of the telescopes in 
x, and whether to use a rectangular or hexagonal 
grid, must still be determined by simulation. 

Us ing the ROOT software package ( Brun and Rademakt .rf 
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19971 ). a Monte Carlo simulation program that in- 
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corporates the effects of atmospheric scintillation 
and camera/ ADU statistics, and allows for the 
variation of placement of telescopes in patterns 
around the satellite ground track, was developed. 
Plots of the average uncertainty, per observation, 
on the measured CALIPSO beam pulse energy 
reaching the ground, for rectangular and hexago- 
nal 4x5 grids of 16" telescopes, as functions of 
the length of the telescope network in the x di- 
rection, are shown in Fig. Q3J It is seen that the 
minimum uncertainty per observation is obtained 
with a hexagonal grid of length 205 m, as shown 
in Fig. [13l As shown in the plot, the average 
uncertainty per observation for this pattern is ap- 
proximately 1.1%. This includes effects of atmo- 
spheric scintillation and camera/ADU statistics, 
but does not include uncertainty due to labora- 
tory and on-site absolute calibration of the tele- 
scopes and cameras, which are effects that at this 
level of uncertainty would be likely to dominate 
in the absence of a very careful program for abso- 
lute calibration. The simulation also includes the 
effect of observations which fail to yield a recon- 
structible measurement of the observed energy for 
a CALIPSO pulse (approximately 4.5% of obser- 
vations with the optimized setup), primarily due 
to the network falling in between (in x), or to the 
left or right (in y), the train of CALIPSO pulses; 
such failed observations contribute zero to the av- 
erage of the inverses of the standard deviations of 
the simulated observations. 

In order to test the operation of a network of 
larger telescopes, a single Meade 16" telescope, as 
shown in Fig. [15J was added to the center of the 
camera and photodiode network. This telescope 
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Fig. 14. — The average uncertainty, per obser- 
vation, on measurements of the CALIPSO beam 
pulse energy reaching the ground (as a fraction 
of that energy) for (top) a rectangular 4x5 grid 
of twenty 16" telescopes and (bottom) a hexagonal 
grid of twenty 16" telescopes (as shown in Fig. [T3J) , 
as functions of the length of the grids in the direc- 
tion of the CALIPSO satellite motion. This uncer- 
tainty includes effects due to atmospheric scintilla- 
tion and camera/ADU statistics, but does not in- 
clude absolute calibration or beam profile effects. 
The values were determined using 5000 MC "ob- 
servations" for each of the twenty bins included on 
each plot. As seen above, the optimum grid pat- 
tern for the 20 telescopes is a hexagonal grid 205 
meters in length. 



has been sucessfully operated, together with the 
network of the seven small Panasonic cameras, at 
a single CALIPSO overpass near Comox, B.C., in 
July of 2009. 
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Fig. 15. — Calibration dark box for the 16" tele- 
scope. A 2-D moveable stage for an attenuated 
532 nm laser beam for calibration is on the top of 
the box. 



7. FUTURE OF THE GROUND-BASED 
OBSERVATION NETWORK 

As shown above, the uncertainties on atmo- 
spheric extinction in the measurements using the 
Panasonic camera network are dominated by at- 
mospheric scintillation, with the second-largest 
source being beam profile-related uncertainty. 
Both of these uncertainties can be dramatically 
reduced by using larger-aperture devices (as de- 
scribed in the MC studies above), and by pre- 
cise pre-flight measurement of beam profile shape. 
However, moving a network of 16" telescopes poses 
some logistical challenges that are not present 
with a network of small consumer cameras. Stor- 
age, transport, and setup facilities are required 
with a number of larger devices. Such challenges 
are certainly not insurmountable, though. As an 
intermediate step, a smaller network of six 16" 
telescopes, together with the present Panasonic 
cameras, is more easily transportable, and would 
reach a precision of 6.0% average uncertainty per 



observation (when the significantly increased level 
of observations that fail to yield a reconstructable 
measurement is taken into account), a little over 
a factor of 5 less precise than a twenty telescope 
network. 

We envision a small storage area for the tele- 
scopes and other hardware in the southwestern 
U.S., with facilities for moving and setup of the 
telescope network below an expected satellite over- 
pass, with each setup for a twenty 16" telescope 
network as shown in Fig. [13] As LIDAR-based 
satellites for atmospheric science will continue 
to be available following the completion of the 
CALIPSO mission (see below), and as the geo- 
location of laser pulses is important for such mis- 
sions (in addition to the motivation of atmospheric 
calibration), it is envisioned that the time and 
manpower necessary for such a facility will become 
available. 

8. FUTURE OF SATELLITE-BASED 
CALIBRATION STANDARDS 

Following the completion of CALIPSO's mis- 
sion, it is envisioned that, as per recommenda- 
tions of the National Academy of Sciences, a 
follow-up mission, known as ACE (Aerosol-Cloud- 
Ecosystems), will succeed CALIPSO, and continue 

to provide a well-cali brated source of visible light 

from low Earth orbit (|NAS Committee on Earth Science 
2008). ACE is forseen to have an increase in pulse 
frequency above that of CALIPSO, to 100 Hz, 
which will improve the ability to measure pulses 
from the ground, as the pulses will not be so widely 
separated. Thus, the ground-based network can 
be made more compact, and/or one will have the 
ability to measure multiple pulses per overpass, 
increasing the data sample. With a network of 
twenty 16" telescopes as shown above, we expect 
to obtain uncertainties on extinction, and thus 
absolute photometry, to better than 1.1% for mul- 
tiple pulses from each observation of ACE, limited 
by only by the precision of laboratory and on-site 
absolute calibration of the telescopes and cameras. 

However, such measurements would still be 
valid only for the zenith, at the specific location 
of the telescope network, for the specific laser fre- 
quency. To obtain results that apply more directly 
to calibration of photometry for the majority of as- 
tronomical results, the satellite laser source should 



21 



be able to point at major observing facilities, and 
modify its frequency to be able to scan through the 
visible and near-infrared spectrum. Both capabil- 
ities are definitely within reach of technology, as 
lasers that are tunable from the near-IR, through 
the visible, through near-UV, as well as mov- 
able mirrors to point a beam, are presently com- 
mon features of laboratories and could be qual- 
ified for space. However, such capabilities will 
require the generation of hardware beyond ACE. 
Both pointability (to scan over temporal atmo- 
spheric features, e.g. hurricanes, forest fires, and 
man-made sources of pollution), and wavelength 
tunability, are of interest to CALIPSO and ACE 
atmospheric scientists as well. 

Beam-based pulsed sources are not the only op- 
tions for future photometric calibration standards 
above the atmosphere, however. More isotropic, 
and continuous, sources have the advantage that 
measurements are far less sensitive to the precise 
relative angle of the source and the observer, and 
that one could time-average over atmospheric scin- 
tillation, rather than having to aperture-average. 
Thus, a more isotropic calibrated source above 
the atmosphere is an attractive option. No such 
sources presently exist. In order to test such equip- 
ment for a future satellite of that type, and to ob- 
tain valuable data in the interim, the author and 
collaborators intend to fly calibrated light sources 
above observatories on a high-altitude balloon. A 
very preliminary sketch of the payload we intend 
to fly may be seen in Fig. [16] Such balloon flights 
would not only provide important data for testing 
future satellite equipment, they would also have 
the advantage over satellite sources that the pay- 
load may be recovered following the flight, and 
tested in the laboratory, rather than only having 
such laboratory tests preceeding the launch. How- 
ever, they of course cannot attain the altitude, nor 
can they easily attain the global reach, of a satel- 
lite. 

One must note, however, that any future 
satellite-based light source would still clearly 
occupy a fixed orbit, and thus a rigorous pro- 
gram of atmospheric monitoring (using techniques 
such as ground-based elastic LIDAR backscat- 
ter, an d dual-band GPS for water vapor mon- 
itoring 
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be required to calibrate the full sky. A device 
that could provide an ability to calibrate any se- 
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Fig. 16. — Conceptual sketch of a high altitude 
balloon-borne lamp calibration system for obser- 
vatories. The concept is based on four integrating 
spheres, each sphere containing either a quartz- 
tungsten-halogen (QTH) lamp, or a light-emitting 
diode (LED) or plasma lamp. An electrical substi- 
tution radiometer (ESR), along with silicon trap 
photodiodes, provides absolute radiometry of the 
sources. 



lected point in the sky at chosen times, as well as 
the features of visibility from major observatories 
and wavelength-tunability, is a far more distant 
concept, althoug h recent progr ess in reusable sub- 
orbital vehicles (|NASAl l2010h could potentially 
point the way toward such an eventual possibility. 
Even then, present and traditional means of atmo- 
spheric calibration would likely remain necessary. 
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In any case, the addition of man-made calibrated 
light sources in space to the arsenal of techniques 
for photometric calibration will, however, provide 
a powerful new tool for increasing precision in 
astrophysics. 

9. CONCLUSION 

In a campaign of observations of the 532 nm 
laser beam from the CALIPSO satellite, using a 
network of small cameras and photodiodes, we 
have obtained measurements of atmospheric ex- 
tinction at four dates and locations. While the 
uncertainties are presently large, they are dom- 
inated by sources that allow dramatic improve- 
ment, to the (9(1%) uncertainty level, with future 
data. Uncertainties due to atmospheric scintilla- 
tion can be greatly reduced by using optics with 
larger aperture. Beam profile-related uncertain- 
tics can be greatly reduced by precise pre-flight 
measurement. A network of multiple 16" tele- 
scopes, together with present and upcoming LI- 
DAR satellites, will allow such improved datasets 
to be obtained. 

Additionally, measurements of extinction using 
a more isotropic source are planned, which avoid 
the requirement of simultaneous precise measure- 
ment of beam location as well as observed power. 
Such equipment will initially be tested on high- 
altitude balloon flights. 

Improved precision in photometric calibration 
will be nearly as critical for astronomy as increased 
aperture telescopes in upcoming decades. The fu- 
ture of precision photometry is extremely promis- 
ing, and laboratory-based standards in space, such 
as described above, allow one to forsee many-fold 
improvement in photometric calibration as a near- 
term prospect. 
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